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Abstract. We discuss the relevant progress that has been made in the last few years 
on the microscopic theory of the pairing correlation in nuclei and the open problems 
that still must be solved in order to reach a satisfactory description and understanding 
of the nuclear pairing. The similarities and differences with the nuclear matter 
case are emphasized and described by few illustrative examples. The comparison of 
calculations of different groups on the same set of nuclei show, besides agreements, 
also discrepancies that remain to be clarified. The role of the many-body correlations, 
like screening, that go beyond the BCS scheme, is still uncertain and requires further 
investigation. 
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1. Introduction 

Despite tlie pairing in nuclei was established more than fifty years ago, the underlying 
interaction processes that are responsible of this remarkable and important correlation 
are not yet understood completely. The main reason of the difficulty is that the 
microscopic theory of nuclear pairing must rely on the effective interaction among 
nucleons, that is not known from first principles. Furthermore it is strong and therefore 
it requires the solution of a complex many-body problem. Even in nuclear matter, where 
the problem is hoped to be simplified to some extent, an accurate microscopic theory is 
not available. Numerical Monte-Carlo "exact" calculations are present in the literature 
only for low-density pure neutron matter [1]. Unfortunately these Monte-Carlo estimates 
are not in satisfactory agreement among each others [1] and in any case they can hardly 
elucidate the microscopic mechanisms which are at the basis of the onset of pairing. 

Besides the many-body aspects of the problem, at least other two features of the 
nuclear pairing have to be mentioned. One is related to the fact that the pairing 
phenomenon occurs close to the Fermi surface, while the bare nucleon-nucleon (NN) 
potential necessarily involves also scattering to high energy (or momentum) due to its 
strong hard core component, which is one of the main characteristics of the nuclear 
interaction. It looks therefore natural to develop a procedure which removes the high 
energy states and "renormalize" the interaction into a region close to the Fermi energy. 
This can be done in different ways, among which the most commonly used seems to 
be the Renormalization Group (RG) Method [2]. A second feature is the relevance of 
the single particle spectrum, not only because the density of states at the Fermi surface 
plays of course a major role, but also because the whole single particle spectrum has 
infiuence on the effective pairing interaction. 

Despite these uncertainties of the microscopic theory of pairing in nuclear matter, 
there is a commonly accepted point of view that the gap value A is quite small at the 
normal density po, much smaller than typical values of heavy atomic nuclei. This is an 
indication that the nuclear surface must play a major role in establishing the value of 
the pairing gap in finite nuclei. 

In the last few years relevant progresse have been made in the microscopic 
calculations of pairing gap in nuclei [31 IH El El [3 El El [in]. The main established results 
is that the bare NN interaction, renormalized by projecting out the high momenta, is a 
reasonable starting point that is able to produce a pairing gap which shows a discrepancy 
with respect to the experimental value not larger than a factor 2. In view of the great 
sensitivity of the gap to the effective interaction this result does not appear obvious. 
The effective pairing interaction constructed within this renormalization scheme pT| 
explains qualitatively also the surface relevance. Indeed, this interaction at the surface 
can exceed the value inside by one order of magnitude. Direct effect of the surface 
enhancement of the gap was presented in jl2[ [13] for semi-infinite nuclear matter and 
in [H] for a nuclear slab. 

The role of the surface is also apparent in explaining the puzzling value of the 
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Figure 1. Cooper pair distribution for the ^^"Sn nucleus. The self-consistent mean 
field is found with the Energy Density Functional DF3 of [TB] ■ 



coherence length extracted from the pairing gap. In fact, for a pairing gap of the order 
of 1 MeV, the coherence length is expected to be much larger than the size of the nucleus. 
It can be argued that at the surface the pairing gap is larger and the coherence length 
shorter, i.e. the size of the Cooper pairs should shrink. This was explicitly shown in a 
recent paper by N. Pillet, N. Sandulescu, and P. Schuck [15] on the basis of the HFB 
approach and the effective DIS Gogny force. It was shown that Cooper pairs in nuclei 
preferentially are located in the surface region, with a small size (2 — 3 fm). In [T6] . 
it was examined to what extent the effect found in [15] is general and independent of 
the specific choice of NN force. This investigation was carried out for a slab of nuclear 
matter, and the pairing characteristics obtained with the Gogny force were compared 
with those with the realistic Paris and Argonne vig forces. The results obtained with 
the two realistic forces agree with each other within an accuracy of about 10% and agree 
qualitatively with those of the Gogny force. In [T^ an analogous analysis, with the Paris 
potential, was made for the nucleus ^^°Sn which is a standard polygon for examining 
nuclear pairing. Again the results turned out to be very close to those of [15]. One of 
them is shown in figure 1, where the "probability" distribution of Cooper pairs, 

p{R) = AttR^ j x\R, r) £r, (1) 

is displayed. Here R, r are cm. and relative coordinates. In Fig. 1 this quantity is 
normalized to the total number of Cooper pairs, A^^cp = / p{R)dR ~ 10. One can see 
that Cooper pairs, indeed, are strongly concentrated on the surface. 

Besides the renormalization of the bare interaction of the high momentum 
components, other physical effects should be included in a microscopic approach, like 
the ones related to the effective mass, or more generally, the single particle spectrum 
and the many-body renormalization of the pairing interaction. The role of these and 
other features of the microscopic scheme are still to be clarified. 

In this paper we will try to summarize the recent achievements and the main open 
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problems that hinder the development of an accurate microscopic many-body theory of 
the nuclear pairing. 

2. Renormalization of the high energy components 

The pairing gap in nuclei is of the order of 1-2 MeV, while the typical energy scale of the 
bare NN interaction is of few hundreds MeV. This indicates the difficulty of controlling 
the microscopic construction of the effective pairing interaction V^g at the Fermi surface. 
In the BCS approach a simplified estimate is commonly used, the so called weak coupling 
limit. In this approximation the relationship between the gap Ap at the Fermi energy 
and the effective interaction has the form 



where is the Fermi energy and up — m*k-p/T:'^ is the density of state at the Fermi 
energy. The exponential dependence makes the problem particularly delicate, since a 
small change in the effective interaction can result in substantial change of the pairing 
gap. In phenomenological theories, like the Finite Fermi System (FFS) theory or the 
Hartrcc-Fock-Bogolyubov method with effective forces, the value of V^cff is considered as a 
phenomenological parameter (or a set of parameters) to be used to fit the data. Despite 
the undoubt success of the phenomenological approaches, the challenge of ah initio 
evaluation of the pairing gap remains one of the most fundamental and not completely 
solved problems in nuclear physics. 

Let us first consider the problem of reducing the interaction to an effective one 
close to the Fermi energy. In general language this can be seen as a typical case 
that can be approached by an "Effective Theory", where the low energy phenomena 
are decoupled from the high energy components. In this procedure the resulting low 
energy effective interaction is expected to be independent of the particular form of the 
high energy components. However the procedure is not unique. The RG method has 
been particularly developed for the reduction of the general NN interaction keeping 
the deuteron properties and NN phase shifts up to the energy where they are well 
established. In this way a potential, phase equivalent to a known realistic NN potential, 
can be obtained, which contains only momenta up to a certain cut-off. The extension of 
this procedure to the many-body problem appears in general to require the introduction 
of strong three-body forces. It has been applied also to the pairing problem. To illustrate 
the difficulty of decoupling low and high momentum components for the pairing problem, 
we consider the simple case of nuclear matter in the BCS approximation. The BCS 
equation for the gap A(A;) can be written as 



where V is the free NN potential, £{k) = e{k) — fi, e{k) is the single particle spectrum 
and /X, the chemical potential. It is possible to project out the momenta larger than a 



A- 



2£Fexp(l/l/Fyeff) , 



(2) 




(3) 
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cutoff kc by introducing tlie interaction Ves{k, k'), wliicli is restricted to momenta k < k^- 
It satisfies the integral equation 

Kfr(fc, k') = Vik, k') - Y: ^^^^^'Ye^^^''^^'^ ' (4) 

k">kc 

where E(k) = -^/£(fcp~+~Ap?j2^. The gap equation, restricted to momenta k < kc and 
with the original interaction V replaced by Ves, is exactly equivalent to the original gap 
equation. The relevance of this equation is that for a not too small cut-off the gap A{k) 
can be neglected in E{k) to a very good approximation and the effective interaction 
depends only on the normal single particle spectrum above the cutoff kc- In the RG 
approach the low momenta interaction Viow-k{k, k') is constructed in such a way to keep 
the half-on-shell two-body T-matrix in free space and therefore the phase shifts. To the 
extent that the pairing gap is determined only by the phase shifts, the two approaches 
are therefore equivalent. 



3. The single particle spectrum and the effective mass 

The pairing gap value is strongly affected by the density of state at the Fermi level, as 
the weak coupling limit (E]) indicates. In turn, the density of state is proportional to 
the effective mass. As it is well known, one has to distinguish between the so-called k- 
mass, which is generated by the momentum dependence of the single particle potential, 
and the e-mass, which is generated by the energy dependence of the single particle self- 
energy. The total effective mass is just the product of these two effective masses. The 
inclusion of an energy dependent single particle self-energy produces also the so-called 
Z-factor, i.e. the quasi-particle strength. 

In nuclear matter the four-dimensional gap equation incorporating the momentum 
and energy dependent irreducible particle-particle vertex V(fc, e; fc', e') and the self- 
energy S(fc, e), reads [T9| [201 [2T] 



with 

D{k, e) = {gg')-^ = [M{k, +e)-e- iO] [M{k, -e) + e - iO] + A^k, e) (6) 

and 

M{k,e) = — + E{k,fi + e)-fi, (7) 
2m 

where we define for convenience the energy e relative to the chemical potential /i. 

and ^ in ([6]) are, correspondingly, the single particle Green functions with and 
without pairing. For realistic systems, V and S cannot be determined in exact way 
and significant approximations have to be performed. The usual BCS approximation 
amounts to replacing the interaction vertex by the (energy independent) bare nucleon- 
nucleon potential V, and the nucleon self-energy by some realistic s.p. spectrum. The 
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latter is characterized by the single particle potential, that can have, in nuclear matter, 
a momentum dependence which persists even at high momenta. To illustrate this point 
we have reported in Fig. |2] the effective mass obtained within the BHF approximation 
as a function of momentum at different Fermi momenta. The effective mass has usually 
a peak close to the Fermi momentum, but it is substantially different from the bare one 
even at high momenta. 

As shown in detail in Ref. [21j, energy dependence of the self-energy can be taken 
into account explicitly provided the interaction vertex remains static, V{k,e; k',e') = 
V{k, k'). In this case, the general gap equation (|6]) is reduced to the form 

Am = - V ^MEELAf A:') (81 

which reminds the BCS gap equation ([3]), with the "symmetrized" s.p. energy 
M,(k) ^ Re ^M(k,+e,)+^M(K-e,)-^ 

appearing in the denominator and with the kernel modified by the spectral factor 

In table [U the neutron effective masses at the Fermi surface m*{kp) of symmetric 
and asymmetric, {N — Z)/A = 1/6, nuclear matter found within the BHF approach 
at the equilibrium density po is compared with those of Skyrme forces. In the second 
line, the difference S{m* /m)s.nm = {fn* /^)a.nm — ijn* I'm) snm is given. As one can see, 
the SLy4 effective mass is rather close to the BHF one at the Fermi surface. Other 
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Figure 3. Effective mass rrip at the Fermi surface depending on the Fermi momentum 
fcp- 

Table 1. Comparison of the neutron effective mass in the symmetric and asymmetric 
nuclear matter from the BHF calculation with that of different kinds of the Skyrme 
force. 





SKP [H] 


SKM* P3' 


SLy4 [M] 


BHF 


{m*/m)snm 


0.85 


0.71 


0.74 


0.78 


(5(TO*/m)a„m 


0.22 


0.12 


-0.04 


-(0.03 0.05) 



two kinds of Skyrme force give absolutely different isotopic asymmetry dependence of 
m*. The density dependence of the effective mass in symmetric nuclear matter for 
BHF and the considered Skyrme forces is reported in figure [31 Around saturation the 
trends are smooth, but at lower density the BHF effective mass has a behavior that 
interpolates between different Skyrme forces, and therefore it cannot be reproduced by 
a given Skyrme force. 

To illustrate the relevance of the high momenta components k > kc on determining 
the effective pairing interaction V^g around the Fermi momentum, we present in figure 
m the value of Ves at the Fermi momentum obtained by solving Eq. (jll). We take 
kc = V2kp and consider three cases : i) the free spectrum, ii) the spectrum obtained 
from the BHF calculation, and iii) the spectrum with a constant effective mass, taken 
at the Fermi momentum and from the BHF results. The calculations are performed in 
the density range relevant for the bulk and surface regions of finite nuclei. 
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Figure 4. Effective pairing interaction for the free single particle spectrum (dotted 
line), the BHF spectrum (full line) and the constant effective mass spectrum (dashed 
line). 





Figure 5. Pairing gap at the Fermi momentum for the three cases of free single 
particle spectrum considered in Fig. (j4|. 



The substantial difference between cases i) and iii) shows the relevance of the high 
momentum components of the single particle spectrum. On the other hand, the strong 
overlap of the curves corresponding to the cases ii) and iii) suggests that the details of 
the momentum dependence of the effective mass are less relevant. This will be helpful 
for the analysis in finite nuclei, where the approximation of a constant effective mass 
will be made. The corresponding pairing gaps at the Fermi momenta are reported in 
figure [51 

4. The many-body problem 

In atomic nuclei, a new branch of low-laying collective excitations appear, the surface 
vibrations ("phonons"). They could be interpreted as a Goldstone mode corresponding 
to spontaneous breaking of translation invariance in nuclei [25]. The ghost dipole 
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phonon is the head of this branch. In the self-consistent FFS theory or other self- 
consistent approaches, say, the SHF+RPA (or SHFB+qRPA) method, the constraint 
w^- = is fulfilled identically and the next members of the branch of natural parity 
states, 2f, 3^ and so on, have small excitation energies, less of a typical distance between 
neighboring shells, uq — €p/A^^^. Up to now, only phenomenological approaches were 
used to describe the surface vibrations and their role in the gap equation via the induced 
interaction. The latter was examined by the Milano group in a series of papers using 
various approaches, with the application to the nucleus ^^°Sn. In Ref. [1] the study 
was performed within the Nuclear Field Theory (see [26] and References therein) in the 
two-phonon approximation. One and two-phonon terms of the induced interaction were 
taken into account, as well as the corresponding corrections to the single-particle states. 
The latter includes both shifts of the single-particle energies and spread of their strength 
which is analogous to the Z-factor in (IHl). The "direct" term of the gap, Adir — 0.7 MeV, 
was found in [4j with the Argonne V14 NN-potential and the mean field generated by the 
SLy4 Skyrme force with the effective mass m* ~ 0.7m. It is one half of the experimental 
gap (in which is estimated as A^xp — 1-4 MeV, while the simplest 3-point formula 
yields Aixp^ ^ 1.3 MeV). By including all the corrections due to low-laying surface 
vibrations, the complete value A = Ajir+ind turns out to be in agreement with the 
experiment. 

In [27] the induced interaction itself was analyzed in detail, again concentrating on 
low-lying {ul < 5 MeV) surface vibrations. The Viow-k potential was used as the free 
NN-interaction, instead of Argonne V14 in [1], and again the experimental value of the 
gap was obtained as a sum of practically equal "direct" and "induced" contributions 
multiplied by the Z-factor representing the quasiparticle strength at the Fermi surface. 
The latter was estimated on the basis of calculations in [1] as Z = 0.7. 

An essentially different approach was used in [5] where all collective states with 
spin values J < 5 and excitation energies uj < 30 MeV were included into the induced 
interaction term. The surface vibrations discussed above are only a part of them. Both 
the natural and unnatural parity states contribute to the sum. The qRPA method with 
the SLy4 Skyrme force was used for every J'^-channel. The fragmentation and self- 
energy effects were approximately taken into account with the following relation for the 
total effective pairing interaction: 

< 12|ydir+ind|34 >=Z< 12|Kiir + ^nd|34 > , (11) 

where Z = 0.7 value again was used. A short notation |12 > is used in (fTTl) for the two- 
particle basis states. To make easier the analysis and comparison with other calculations, 
separate calculations of the gap from Vdir (Argonne V14) and Vi^d, both without the Z- 
factor, were made in [5]. They give Adir = 1.04 MeV and Aind = 1.11 MeV. Using the 
recipe of Eq. (fTTj) . yields the average gap Ap = 1.47 MeV which is again close to the 
experiment. 

In our opinion, each method used in these papers has some weak points which could 
be criticized. In the first case, the contribution of so-called tadpole diagrams [25[ [28] 
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was neglected, see also discussion in |21]. In the last case, the use of Skyrme force in the 
particle-hole spin and spin-isospin channels (unnatural parity states) is questionable. 
Indeed, the Skyrme parameters were fitted only to phenomena which are related to 
"scalar" Landau-Migdal amplitudes, F,F'. The combinations which determine the 
spin-dependent amplitudes G, G' were not checked up to now by an attempt to describe 
corresponding phenomena, e.g. magnetic moments, Ml-transitions, and so on. This 
problem is discussed also in [5], and some complementary calculation was made with 
G = G" = 0. It resulted in a very strong induced interaction yielding very big gap 
Ap = 2.12 MeV. This could be considered as an estimate of the uncertainty of such 
calculations with the use of phenomenological forces. Indeed, even the amplitudes 
F, F' which are known sufficiently well in vicinity of the Fermi surface could change 
significantly when high energy excitations are considered. Another questionable point 
is the use in this calculation scheme of the Z-factor found only from the low-lying surface 
vibrations. The high energy response function included into the V^nj will contribute to 
the Z-factor as well. This contribution comes mainly from the spin-isospin channel and 
could be estimated as ^ 0.8 [30] . 

We see that the problem of the screening effect is, indeed, very difficult, and some 
approximations must be introduced in the calculations. In our opinion, the fact that 
essentially different methods were used in [1] and [5] with different results for Aind 
shows by itself that the problem of finding contribution of the induced interaction into 
the pairing gap in atomic nuclei is far from being solved. 

5. Solution of the "afe initio'''' gap equation in finite nuclei 

The explicit form of the gap equation ([5]) in the coordinate representation for a non- 
uniform system is as follows [31] : 



As in Sections 2-4, single-particle energies e, e' are counted off the chemical potential 
/i. Dealing with nuclear matter, we set yU = yUg — ~16 MeV (the leading term in the 
Weizsaecker mass formula), whereas we have /i ~ —8 MeV for stable atomic nuclei, as 
the ^^°Sn nucleus considered below. As it was discussed above, in this Section we limit 
ourselves with the simplest Brueckner-like approach in which the irreducible vertex V 
coincides with the free A^A^-potential, V = V , which is independent of energy. In 
this case, the gap A is also independent of energy; hence, the product of two Green's 
functions in ( !T2|) can be integrated with respect to e': 



A(ri, ra, ^) = j "l^l^i, ^a, rg, E = 2/i, e, e') x 

de' 

X ^(rg, rs, e')^'(r4, re, -e')A{r5, rg, e')—dY^dYidY^dYQ. 



(12) 




(13) 



The gap equation (fT2l) can be written in a compact form as 



A = VA'A. 



(14) 
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Below we deal with the BCS gap equation, not the general one, Eq. (IT^ . This explains 
why the term ab initio in the title of the Section is in quotes. In fact, we speak just 
about the first step into the problem, i.e. the solution of the BCS-like gap equation 
with the free A^A^-potential as the pairing interaction. The mean field potential (or 
more general, the mass operator) used in this solution is taken as a phenomenological 
input. To go beyond the BCS approximation it is necessary, within an ab initio method, 
to calculate, first, the mass operator and, second, corrections to the interaction block V 
in the gap equation ( |T2|) . The latter includes the induced interaction discussed above 
(see also Ref. [32]) and three-body forces [33]. It turns out that, even at the level of 
this simplest ab initio calculations, serious contradictions are still present. 

The integral equation f|T^ can be reduced to the form adopted in the Bogoliubov 
method. 



where the abnormal density matrix x = A^A can be expressed in terms of u , f-functions. 



which satisfy the system of Bogoliubov equations. The summation in (fT6ll is over the 
complete set of Bogoliubov functions with the eigen energies Ei > 0. 

The Milano group was the first who, in a series of papers |3|, Hj, [5] and Refs. 
therein, solved the gap equation ( |T5l) with the realistic Argonne NN-force Vi4 for the 
nucleus ^^°Sn. The latter was chosen for a definite reason. Indeed, the chain of semi- 
magic tin isotopes is a traditional polygon for examining nuclear pairing [22], |18j . 
The nucleus under discussion is in the middle of the chain and the number of neutrons 
participating in the pairing, those above the closed shell = 50, is sufficiently big to use 
the approximation of the "developed" pairing, used, in fact, in ( |T5i) . This approximation 
implies neglecting the particle number fluctuations [31] typical of the BCS-like theories. 
The set of Bogoliubov equations was solved directly in the basis {A} of the states with a 
fixed limiting energy -Emax- Such direct method is difficult because of a slow convergence 
of sums over intermediate states A in the gap equation. These sums are analogous to 
integrals in the momentum space in the gap equation for infinite nuclear matter, see 
Section 2. In [3] the value of E^^^=QOO MeV was used, and in [U [5], E^^^=800 MeV. 
The analysis of [7] showed that the use of such big value of -Emax permits to find A only 
with accuracy of 10%. In [3j the Shell Model basis was used with the Saxon- Woods 
potential and the bare mass, m* = m, and the value A=2.2 MeV was obtained which is 
by a factor one and half greater than the experimental one (~1.3-i-1.4 MeV). Evidently, 
this contradiction forced the authors to use in further works the self-consistent basis 
of the Skyrme-Hartree-Fock (SHE) method with the density depending effective mass 
m*{p) 7^ m. In particular, the popular Sly4 force was used, which is characterized by 
a small effective mass, equal to m* ~ 0.7m in nuclear matter at the normal nuclear 
density. Solving the BCS gap equation with such a basis the value of A ~ 0.7 MeV was 
obtained in [1] and ~ 1.0 MeV in [5]. A close value for the gap in the BCS equation 



A = -Vx, 



(15) 




(16) 
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was found in [7j for the nuclear slab with parameters which mimic ^^'^Sn nucleus. This 
calculation was done with Argonne vig potential which differs only slightly from the 
Vi4 one, but the single-particle basis with m* = m was used. Remind that in [H |5] 
corrections to the BCS gap due to induced interaction were considered bringing the 
results rather close to the experimental data as discussed in Section 4. The crucial 
dependence of the gap on m*, which is easily seen in weak coupling limit (|2]) for the 
gap in nuclear matter, is, of course, the main reason of so strong variation of the BCS 
gap from [3] to [5]. In finite nuclei, this dependence is weaker than in nuclear matter, 
as the surface region plays here the main role and one has m*(p) — )■ m in this region. 
However, the m* effect remains strong. 

Recently, results from the ab initio BCS equation f|T5l) for a number of semi-magic 
nuclei were published [9l[10]. They are based on the soft realistic low-k force discussed 
above which was calculated starting from the Argonne vig potential, with the same self- 
consistent Sly4 basis, i.e. the same effective mass, as in [5]. For the nucleus ^^°Sn under 
consideration the value A~1.6 MeV was obtained. This value exceeds the experimental 
one, which raises some questions. Indeed, although there are discrepancies in absolute 
value of the corrections to the BCS gap (see for instance [5] and [35]), their sign is 
more or less definite. All calculations of these corrections, to our knowledge, increase 
A considerably. Hence the BCS equation has to lead to the gap value smaller than the 
experimental one. In addition, there is a direct contradiction between the results of the 
Milano group and the ones of Duguet with co-authors, despite the BCS problem was 
solved with very similar inputs. Evidently, there is some difference in the method of 
including the effective mass, which is hidden in [9], [10] under the renormalization of Vig 
into Viow-k- 

In [16] an attempt was taken to clarify the reasons of this contradiction. The BCS 
gap equation f|T5l) was solved for the same ^^*^Sn nucleus with the separable form of 
Paris potential, which simplifies calculations. Our experience of calculations for nuclear 
slab [HI [7] shows that the difference between the Paris potential and the Argonne vig 
one for the gap value is of the order of 0.1 MeV which is significantly smaller than 
the deviation under discussion. For projecting out the high momenta contribution, the 
so-called Local Potential Approximation (LPA) method was used. This new version of 
the local approximation was introduced by our group for semi-infinite nuclear matter 
and nuclear slab system, see the reviews [29l [36] . In general, this method is analogous 
to the renormalization scheme for solving the gap equation in infinite nuclear matter 
described in Section 2. To solve the gap equation in the form f[T4l) for finite systems, 
we split the complete Hilbert space S of the pairing problem into the model subspace 
5*0, including the single-particle states with energies less than a fixed value Eq, and the 
subsidiary one, S'. Correspondingly, the two-particle propagator ([T3!l is the sum 

A' = A'q + A'. (17) 
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The notation becomes obvious if one expands fll3p in the basis of single-particle functions 
0A(r), 

A^(ri,r2,r3,r4)= J] ^^^,^,3,^ 01^(ri)0I,(r2)0A3(r3)K(r4) • (18) 
A1A2A3A4 

The model space propagator Aq includes the terms of the sum f fTSj) with single-particle 
energies ex < Eq, A' being the remaining part. 

The gap equation is reduced to the one in the model space, 

A = KffAgA, (19) 

with the effective pairing interaction obeying the integral equation in the subsidiary 
space, 

VeS=V + VA'VeS. (20) 

The LPA method concerns the solution of (!20l) . Solving this equation directly in 
coordinate space is rather complicated, not much simpler than the gap equation (fT^ 
in the complete Hilbert space. The problem could be simplified with the use of the 
LPA. It turned out that, with very high accuracy, for each point R, one can use the 
formulas for the infinite system in the potential field f/(R) (it explains the term LPA). 
This simplifies equation for V^s significantly, in comparison with the initial equation for 
A. As a consequence, the subspace S' could be chosen as large as necessary. Validity 
of the LPA could be justified by finding that, beginning from some value of Eq, the 
result for the gap doesn't change under additional increase of Eq. It turns out that 
for a nuclear slab, the value of £"0— 20-=-30 MeV is sufficient P, [7j. For finite nuclei, it 
should be taken a little bigger, Eq = 40 MeV [8]. In the slab calculations, we used the 
bare nucleon mass, m* = m.ln this case, the reduction of high momenta components of 
the NN-force in fl20|) is, in fact, quite similar to the renormalization procedure resulting 
in a low-k interaction [37]. The difference is that in the LPA the renormalization is 
coordinate dependent. 

The gap equation in the model space 5*0 was solved in the A-representation with 
the use of different single-particle bases (px. A discretization method of the continuum 
spectrum was used with the wall radius L=16 fm. Increasing the radius to L=24 fm does 
not practically influence the results. The radial eigen-functions Rnij{r) were found with 
a step h=0.05 fm. We have used the Shell-Model basis with the Saxon-Woods potential 
with a standard set of parameters and also with several self-consistent basis obtained 
with different methods : the Generalized Energy Density Functional method by Fayans 
et al. [18] with the functional DF3, and the SHF method as well with different kinds of 
Skyrme forces. The bare mass, m*=m, is used in the first method, just as in the Shell 
Model, whereas in the SHF method the effective mass is not equal to m and is density 
dependent. To clarify the role of the effective mass, we chose two kinds of the Skyrme 
force, SKP and SKM*, for which the difference between m* and m is quite moderate, 
and the popular SLy4 force, with significant difference of m* from m, see Fig. [61 We 
see that the SKP effective mass deviates significantly from that for symmetric nuclear 




Figure 6. Coordinate dependence of tlie effective mass for different Skyrme forces in 
the ^^"Sn nucleus. 



matter, Fig. [21 This is due to the strong isospin asymmetry effect {{N — Z)/A = 1/6 for 
^^°Sn nucleus) for this kind of Skyrme force. For SLy4 and SKM* this effect is rather 
modest. Remind that the Sly4 basis was used in calculations of A by Milano group and 
by Duguet with coauthors as well. As to the calculation of the effective interaction in 
the subspace S", we first put m*=m. In table 2 the diagonal matrix elements A^a of 
the neutron gap in the ^^°Sn nucleus, for 5 levels nearby the Fermi level, are given for 
each basis under consideration. The quantity Ap is the corresponding Fermi-average 
value: Ap = (^J + 1)^aa/ SaI^J + !)• As we see, in all cases except of the last one 
the found value of the gap exceeds the experimental one, 1.4 MeV, significantly. This 
indicates the necessity to take into account the difference between the effective and bare 
masses in the ab initio BCS gap equation. 

In Fig. [71 for each kind of mean field, the anomalous density z^(-R) = x{R, r = 0) 
is displayed. Here the notation R = (ri + r2)/2, r = r2 — ri is used. As the gap A is 
proportional to this quantity, i^{R) depends on the effective mass in the same way as 
the matrix elements Aaa in table 2. Indeed, the anomalous densities for the effective 
forces DF3, SKP and SKM* are rather close to each other, whereas for the SLy4 force 
with a small effective mass the anomalous density is suppressed significantly. There 
is a pronounced surface maximum of i^{R) for all the cases. This leads to the surface 
enhancement of the pairing gap found in [H] for the slab and in [15], for spherical nuclei. 
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Table 2. Diagonal matrix elements Aw (MeV) with the Paris potential for several 
kinds of the self-consistent basis. 



A 


sw 


DF3 


SKP 


SKM* 


SLy4 


3si/2 


1.52 


1.64 


1.55 


1.55 


1.17 


2(^5/2 


1.60 


1.73 


1.68 


1.64 


1.24 


2^3/2 


1.64 


1.80 


1.71 


1.68 


1.26 


I57/2 


1.85 


2.11 


2.02 


1.91 


1.37 


2^11/2 


1.58 


1.79 


1.69 


1.64 


1.18 


Af 


1.65 


1.85 


1.76 


1.71 


1.25 



0.04 



0.03 - 



DF3 

SKP 

SKM* 

SLy4 




10 



Figure 7. Anomalous density for the ^^°Sn nucleus calculated with different self- 
consistent mean fields. 



Until now, in the case of the Skyrme forces, we took into account the difference 
between the effective and bare masses only inside the model space, whereas we put 
m*=m for calculating Ves- At this point, there is a principal difference between our LPA 
method and calculations of [H |5] and [H |T0] , where the SLy4 effective mass was used for 
all the A-states. For a closer comparison, we made a modification of the LPA method, 
which permits to take into account the density dependent effective mass m*(p„,pp) for 
a part of the space S", including momenta k < A, where A is a parameter. Following 
the idea of LPA, with the potentials Un{R), Up{R) , it is natural to determine at each 
point R the densities of each type t = n,p of nucleons with quasi-classical formulas: 
Pr{R) = (A;J(i?))V37r2, k^R) = [2m;(p„(i?), pp(i?))(/i, - Ur{R))Y/\ where /i„,/ip are 
the chemical potentials of neutrons and protons in the nucleus under consideration. For 
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Table 3. Diagonal matrix elements A\\ with the Paris potential for the Sly4 basis 
depending on the way of account for the effective mass in equation for V^g. 



A 


SLy4 


Sly4-1 


Sly4-2 


Sly4-3 


3Si/2 


1.17 


1.07 


0.88 


0.76 


2^5/2 


1.24 


1.13 


0.93 


0.80 


2c^3/2 


1.26 


1.15 


0.95 


0.83 


I57/2 


1.37 


1.23 


0.99 


0.85 


2/I11/2 


1.18 


1.08 


0.88 


0.75 


Af 


1.25 


1.14 


0.92 


0.80 



Table 4. The same as in table |3l but for the Argonne vig potential. 



A 


SLy4 


Sly4-1 


Sly4-2 


Sly4-3 


3Si/2 


1.23 


1.10 


0.83 


0.56 


2^5/2 


1.32 


1.18 


0.89 


0.61 


2ci3/2 


1.34 


1.20 


0.92 


0.63 


I57/2 


1.48 


1.31 


0.96 


0.64 


2/I11/2 


1.27 


1.13 


0.85 


0.57 


Af 


1.34 


1.19 


0.89 


0.60 



the functional SLy4 we made several alternative calculations with different values of A. 
They are denoted as SLy4-l (A=3 fm'^), SLy4-2 (A=4 fm-^) and SLy4-3 (A=6.2 fm"^). 
The first two versions mimic calculations of [HI HO], the latter, of [U E]. The obtained 
gap values are given in table 3. 

Let us now repeat calculations with the SLy4 basis, but for Argonne force vig. 
Results are given in table 4. Comparison with table 3 shows that the difference between 
gap values for the Paris and Argonne force is of the same order (~ 0.1 MeV) as for 
the slab calculations p'J, with the exception of the SLy4-3 version with the cutoff for 
m* in the equation for Ves equal to A3 = 6.2 fm~^. In this case, the difference is 
~ 0.2 MeV. It is worth of noticing that the sign of the difference changes depending 
on A. Namely the Argonne gap exceeds the Paris one for SLy4 (A = 0) and SLy4-l 
(Ai = 3 fm~^) versions, but becomes smaller for SLy4-2 (A2 = 4 fm~^) and SLy4-3 
runs. Such behavior is qualitatively clear. Indeed, the Paris potential is much harder of 
the Argonne one, therefore the relative contribution to V^s of the momentum region, 
say, between A2 and A3, is less than for the Argonne potential. Therefore, for the gap 
equation itself, the role of the corresponding suppression of V^s due to putting m* < m 
is less in the Paris case than in the Argonne one. 

It is rather difficult to make a direct comparison of table 4 with results of [9j. In 
the first column, the effective mass m* ^ m is introduced only in the model space, 
for e\ < Eq = 40 MeV. In the free space, outside the nucleus, it corresponds to the 
momentum cutoff A = 1.4 fm~^; inside the nucleus we have A ~ 2 fm~^.We see that the 
"average" value of A is less a bit of A = 2 fm~^ in [9J. In the second column (SLy4-l) we 
deal with A = 3 fm^^. Thus, we should attribute to the gap value of [9] (A ~ 1.6 MeV) 
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A 


SKM* 


SKM*-1 


SKM*-2 


SKM*-3 


3si/2 


1.58 


1.52 


1.40 


1.29 


2^5/2 


1.71 


1.63 


1.49 


1.38 


2^3/2 


1.75 


1.68 


1.55 


1.43 


lS7/2 


2.01 


1.92 


1.76 


1.62 


2^11/2 


1.72 


1.64 


1.51 


1.40 


Af 


1.78 


1.71 


1.57 


1.44 



an average value of these two columns, A ~ 1.25 MeV, which is noticeably smaller. 
For comparison with [H [5] we should be guided by the last column (SLy4-3), as the 
corresponding value A = 6.2 fm"^ was chosen to reproduce ii^max = 800 MeV from these 
calculations. Again, this correspondence is not literal as it takes place only outside the 
nucleus where we have m* = m. Inside, due to m* ^ m, the same value of -Emax should 
correspond to smaller value of A ~ 5.5 fm~^. Again we should take a value between 
those of the 3-rd column and the last one, closer to the latter. In any case, it will be 
closer to the result of [1] (~ 0.7 MeV) than of [5] (~ 1 MeV). 

The main observation, common to table 3 (Paris) and table 4 (Argonne), is a 
drastic dependence of the gap on the m*{k) behavior in the subsidiary space. In fact, a 
set of calculations with different cutoff A imitates, very roughly, the fc-dependence of the 
effective mass. In the case that the asymptotic limit m*{k) — )■ m occurs sufficiently far 
away, the pairing gap in nuclei does depend on the effect of m* 7^ m at high momenta. 
It is absolutely lost in calculations of A with low-k force found for small cutoff values. 
Evidently, this is the main reason why gap values found in pj by solving the BCS gap 
equation are so high, keeping in mind corrections due to the induced interaction. 

To confirm the leading role of the effective mass in the problem, we made a series 
of analogous calculations for the SKM* force for which the effective mass in the nucleus 
under consideration is much closer to m, than for the SLy4 one, see Fig. [61 The results 
are given in table [5] with the notation similar to that of table HI We see that even for the 
SKM*-3 version with A = 6.2 fm~^ the result does not leave any room for contribution 
of the induced interaction. 

6. Discussion and conclusions 

In this paper we briefly reviewed the status of the microscopic theory of nuclear pairing. 
Although to date there is no consistent theory of nuclear matter pairing, progress has 
been made in developing approximated methods yielding comparable results on the 
density dependence of the pairing gap. However, the results cannot be applied to finite 
nuclei directly because the nuclear surface plays a leading role in nuclear pairing and the 
standard LDA fails at the nuclear surface. As to finite nuclei, some progress also exists 
especially in solving the simplest, in fact BCS-type, " ab initio" gap equation, where 
the free NN potential is used as the effective pairing interaction. Even this equation 
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is not completely microscopic as far as the phenomenological mean field is used. The 
inherent technical problems were overcome, first, by the Milan group |3l HI |5] and more 
recently by other teams [HJ [10] and [I6]. In the first and the last cases, the nucleus 
^^°Sn was considered as "testing sample", whereas in [9] several isotopic and isotonic 
chains were considered. The Argonne Vi4 force was used in [31 HJ |5], the low-k force 
with the cutoff A = 2 fm~^ in |9], and the Paris potential, in [16]. In this paper 
we performed calculations analogous to [16] for the Argonne vig potential which differ 
only a bit from the V14 one. All the results differ from the experimental gap in ^^°Sn, 
Acxp — 1.3 ^ 1.4 MeV, not more than by a factor two which shows relevance of the 
ab initio BCS gap equation as a starting point for the microscopic theory of nuclear 
pairing. However, more detailed comparison shows that there are contradictions even 
at this "first level" of the problem. Indeed, the BCS gap is ~ 1 MeV in [5] and is 
~ 1.6 MeV in |9j, whereas inputs look quite similar. Namely, both the calculations use 
the SLy4 Skyrme force with the effective mass m* ~ 0.7m; the Argonne V14 force is 
used in [5] and the low-k force in [9], but the latter could be obtained from the first 
one with the RGM procedure. The only important difference of the inputs is the size 
of the momentum space where the effective mass contributes: k < A;max — 6 fm~^ in 
[5] and only k < A = 2 fm~^ in [9]. Indeed, the RGM equation is defined for the free 
NN scattering where the equality m* = m is postulated. In fact, we deal with different 
/c-dependence of the effective mass. The equality m* ~ 0.7m takes place for all momenta 
in and only for /c < A, in [9]. This reason of the contradictory results of [5] and [9] 
was discussed in [lO] and was confirmed with the analysis of [16] for the Paris force and 
in this paper, for the Argonne Vig force. We use the so-called LPA method developed by 
us previously in which high momenta are excluded via some renormalization procedure 
which recall the RGM one but is coordinate dependent and permits to introduce in 
the high momentum space the effective mass into the equation for the effective pairing 
interaction. For the same SLy4 basis, we changed "by hands" the size of the space 
where the equality m* ~ 0.7m takes place, putting m* = m outside. Changing this 
dividing point we evolve from the situation close to that of [9] to the one of 0, although 
the correspondence, of course, is not literal. In the first limit, we obtained the value 
of Ap — 1.25 MeV which is noticeably smaller than that in [9]. In the second one, 
we obtained the result closer to that of [4j than of [5]. But, in our opinion, fixing 
some numerical contradictions is not of primary importance. This disagreement can be 
resolved. Much more important is the very high sensitivity of the pairing gap to the k- 
dependence of the effective mass found in our analysis. This dependence could be hardly 
guessed phenomenologically by a lucky Skyrme-type ansatz. Indeed, any such ansatz 
deals with a number, not a function! Therefore an additional microscopic ingredient, 
namely a theory of the fc-dependent effective mass, is necessary even at this first level 
of the problem. 

The next open problem is the inclusion of the many-body corrections to the BCS gap 
equation. The induced pairing interaction due to exchange by virtual surface vibrations 
and other particle- hole excitations is the main of them [U [5] . Up to now, this problem 



Microscopic evaluation of the pairing gap 



19 



was studied only within phenomenological approaches and is far from being solved. Any 
attempt to attack it from first principles hits the absence of the general microscopic 
nuclear theory. 
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